Enhancing the Small-Scale Screenable Biological Space beyond Known Chemogenomics Libraries with Gray Chemical Matter─Compounds with Novel Mechanisms from High-Throughput Screening Profiles.

Phenotypic assays have become an established approach to drug discovery. Greater disease relevance is often achieved through cellular models with increased complexity and more detailed readouts, such as gene expression or advanced imaging. However, the intricate nature and cost of these assays impose limitations on their screening capacity, often restricting screens to well-characterized small compound sets such as chemogenomics libraries. Here, we outline a cheminformatics approach to identify a small set of compounds with likely novel mechanisms of action (MoAs), expanding the MoA search space for throughput limited phenotypic assays. Our approach is based on mining existing large-scale, phenotypic high-throughput screening (HTS) data. It enables the identification of chemotypes that exhibit selectivity across multiple cell-based assays, which are characterized by persistent and broad structure activity relationships (SAR). We validate the effectiveness of our approach in broad cellular profiling assays (Cell Painting, DRUG-seq, and Promotor Signature Profiling) and chemical proteomics experiments. These experiments revealed that the compounds behave similarly to known chemogenetic libraries, but with a notable bias toward novel protein targets. To foster collaboration and advance research in this area, we have curated a public set of such compounds based on the PubChem BioAssay dataset and made it available for use by the scientific community.

Molecular characterization of chronic cutaneous wounds reveals subregion- and wound type-specific differential gene expression.

A limited understanding of the pathology underlying chronic wounds has hindered the development of effective diagnostic markers and pharmaceutical interventions. This study aimed to elucidate the molecular composition of various common chronic ulcer types to facilitate drug discovery strategies. We conducted a comprehensive analysis of leg ulcers (LUs), encompassing venous and arterial ulcers, foot ulcers (FUs), pressure ulcers (PUs), and compared them with surgical wound healing complications (WHCs). To explore the pathophysiological mechanisms and identify similarities or differences within wounds, we dissected wounds into distinct subregions, including the wound bed, border, and peri-wound areas, and compared them against intact skin. By correlating histopathology, RNA sequencing (RNA-Seq), and immunohistochemistry (IHC), we identified unique genes, pathways, and cell type abundance patterns in each wound type and subregion. These correlations aim to aid clinicians in selecting targeted treatment options and informing the design of future preclinical and clinical studies in wound healing. Notably, specific genes, such as PITX1 and UPP1, exhibited exclusive upregulation in LUs and FUs, potentially offering significant benefits to specialists in limb preservation and clinical treatment decisions. In contrast, comparisons between different wound subregions, regardless of wound type, revealed distinct expression profiles. The pleiotropic chemokine-like ligand GPR15L (C10orf99) and transmembrane serine proteases TMPRSS11A/D were significantly upregulated in wound border subregions. Interestingly, WHCs exhibited a nearly identical transcriptome to PUs, indicating clinical relevance. Histological examination revealed blood vessel occlusions with impaired angiogenesis in chronic wounds, alongside elevated expression of genes and immunoreactive markers related to blood vessel and lymphatic epithelial cells in wound bed subregions. Additionally, inflammatory and epithelial markers indicated heightened inflammatory responses in wound bed and border subregions and reduced wound bed epithelialization. In summary, chronic wounds from diverse anatomical sites share common aspects of wound pathophysiology but also exhibit distinct molecular differences. These unique molecular characteristics present promising opportunities for drug discovery and treatment, particularly for patients suffering from chronic wounds. The identified diagnostic markers hold the potential to enhance preclinical and clinical trials in the field of wound healing.

Sustained Trem2 stabilization accelerates microglia heterogeneity and Aß pathology in a mouse model of Alzheimer's disease.

TREM2 is a transmembrane protein expressed exclusively in microglia in the brain that regulates inflammatory responses to pathological conditions. Proteolytic cleavage of membrane TREM2 affects microglial function and is associated with Alzheimer's disease, but the consequence of reduced TREM2 proteolytic cleavage has not been determined. Here, we generate a transgenic mouse model of reduced Trem2 shedding (Trem2-Ile-Pro-Asp [IPD]) through amino-acid substitution of an ADAM-protease recognition site. We show that Trem2-IPD mice display increased Trem2 cell-surface-receptor load, survival, and function in myeloid cells. Using single-cell transcriptomic profiling of mouse cortex, we show that sustained Trem2 stabilization induces a shift of fate in microglial maturation and accelerates microglial responses to Aß pathology in a mouse model of Alzheimer's disease. Our data indicate that reduction of Trem2 proteolytic cleavage aggravates neuroinflammation during the course of Alzheimer's disease pathology, suggesting that TREM2 shedding is a critical regulator of microglial activity in pathological states.

DRUG-seq Provides Unbiased Biological Activity Readouts for Neuroscience Drug Discovery.

Unbiased transcriptomic RNA-seq data has provided deep insights into biological processes. However, its impact in drug discovery has been narrow given high costs and low throughput. Proof-of-concept studies with Digital RNA with pertUrbation of Genes (DRUG)-seq demonstrated the potential to address this gap. We extended the DRUG-seq platform by subjecting it to rigorous testing and by adding an open-source analysis pipeline. The results demonstrate high reproducibility and ability to resolve the mechanism(s) of action for a diverse set of compounds. Furthermore, we demonstrate how this data can be incorporated into a drug discovery project aiming to develop therapeutics for schizophrenia using human stem cell-derived neurons. We identified both an on-target activation signature, induced by a set of chemically distinct positive allosteric modulators of the N-methyl-d-aspartate (NMDA) receptor, and independent off-target effects. Overall, the protocol and open-source analysis pipeline are a step toward industrializing RNA-seq for high-complexity transcriptomics studies performed at a saturating scale.

Single-cell and bulk transcriptomics of the liver reveals potential targets of NASH with fibrosis.

Fibrosis is characterized by the excessive production of collagen and other extracellular matrix (ECM) components and represents a leading cause of morbidity and mortality worldwide. Previous studies of nonalcoholic steatohepatitis (NASH) with fibrosis were largely restricted to bulk transcriptome profiles. Thus, our understanding of this disease is limited by an incomplete characterization of liver cell types in general and hepatic stellate cells (HSCs) in particular, given that activated HSCs are the major hepatic fibrogenic cell population. To help fill this gap, we profiled 17,810 non-parenchymal cells derived from six healthy human livers. In conjunction with public single-cell data of fibrotic/cirrhotic human livers, these profiles enable the identification of potential intercellular signaling axes (e.g., ITGAV-LAMC1, TNFRSF11B-VWF and NOTCH2-DLL4) and master regulators (e.g., RUNX1 and CREB3L1) responsible for the activation of HSCs during fibrogenesis. Bulk RNA-seq data of NASH patient livers and rodent models for liver fibrosis of diverse etiologies allowed us to evaluate the translatability of candidate therapeutic targets for NASH-related fibrosis. We identified 61 liver fibrosis-associated genes (e.g., AEBP1, PRRX1 and LARP6) that may serve as a repertoire of translatable drug target candidates. Consistent with the above regulon results, gene regulatory network analysis allowed the identification of CREB3L1 as a master regulator of many of the 61 genes. Together, this study highlights potential cell-cell interactions and master regulators that underlie HSC activation and reveals genes that may represent prospective hallmark signatures for liver fibrosis.

CellSIUS provides sensitive and specific detection of rare cell populations from complex single-cell RNA-seq data.

We develop CellSIUS (Cell Subtype Identification from Upregulated gene Sets) to fill a methodology gap for rare cell population identification for scRNA-seq data. CellSIUS outperforms existing algorithms for specificity and selectivity for rare cell types and their transcriptomic signature identification in synthetic and complex biological data. Characterization of a human pluripotent cell differentiation protocol recapitulating deep-layer corticogenesis using CellSIUS reveals unrecognized complexity in human stem cell-derived cellular populations. CellSIUS enables identification of novel rare cell populations and their signature genes providing the means to study those populations in vitro in light of their role in health and disease.

YAP, but Not RSPO-LGR4/5, Signaling in Biliary Epithelial Cells Promotes a Ductular Reaction in Response to Liver Injury.

Biliary epithelial cells (BECs) form bile ducts in the liver and are facultative liver stem cells that establish a ductular reaction (DR) to support liver regeneration following injury. Liver damage induces periportal LGR5+ putative liver stem cells that can form BEC-like organoids, suggesting that RSPO-LGR4/5-mediated WNT/ß-catenin activity is important for a DR. We addressed the roles of this and other signaling pathways in a DR by performing a focused CRISPR-based loss-of-function screen in BEC-like organoids, followed by in vivo validation and single-cell RNA sequencing. We found that BECs lack and do not require LGR4/5-mediated WNT/ß-catenin signaling during a DR, whereas YAP and mTORC1 signaling are required for this process. Upregulation of AXIN2 and LGR5 is required in hepatocytes to enable their regenerative capacity in response to injury. Together, these data highlight heterogeneity within the BEC pool, delineate signaling pathways involved in a DR, and clarify the identity and roles of injury-induced periportal LGR5+ cells.

DRUG-seq for miniaturized high-throughput transcriptome profiling in drug discovery.

Here we report Digital RNA with pertUrbation of Genes (DRUG-seq), a high-throughput platform for drug discovery. Pharmaceutical discovery relies on high-throughput screening, yet current platforms have limited readouts. RNA-seq is a powerful tool to investigate drug effects using transcriptome changes as a proxy, yet standard library construction is costly. DRUG-seq captures transcriptional changes detected in standard RNA-seq at 1/100th the cost. In proof-of-concept experiments profiling 433 compounds across 8 doses, transcription profiles generated from DRUG-seq successfully grouped compounds into functional clusters by mechanism of actions (MoAs) based on their intended targets. Perturbation differences reflected in transcriptome changes were detected for compounds engaging the same target, demonstrating the value of using DRUG-seq for understanding on and off-target activities. We demonstrate DRUG-seq captures common mechanisms, as well as differences between compound treatment and CRISPR on the same target. DRUG-seq provides a powerful tool for comprehensive transcriptome readout in a high-throughput screening environment.

Origin of the spectral shifts among the early intermediates of the rhodopsin photocycle.

A combined strategy based on the computation of absorption energies, using the ZINDO/S semiempirical method, for a statistically relevant number of thermally sampled configurations extracted from QM/MM trajectories is used to establish a one-to-one correspondence between the structures of the different early intermediates (dark, batho, BSI, lumi) involved in the initial steps of the rhodopsin photoactivation mechanism and their optical spectra. A systematic analysis of the results based on a correlation-based feature selection algorithm shows that the origin of the color shifts among these intermediates can be mainly ascribed to alterations in intrinsic properties of the chromophore structure, which are tuned by several residues located in the protein binding pocket. In addition to the expected electrostatic and dipolar effects caused by the charged residues (Glu113, Glu181) and to strong hydrogen bonding with Glu113, other interactions such as π-stacking with Ala117 and Thr118 backbone atoms, van der Waals contacts with Gly114 and Ala292, and CH/π weak interactions with Tyr268, Ala117, Thr118, and Ser186 side chains are found to make non-negligible contributions to the modulation of the color tuning among the different rhodopsin photointermediates.

Wagging the Tail: Essential Role of Substrate Flexibility in FAAH Catalysis.

The serine hydrolase, fatty acid amide hydrolase (FAAH), is responsible for the intracellular degradation of anandamide and other bioactive fatty acid ethanolamides involved in the regulation of pain, inflammation, and other pathophysiological processes. The catalytic site of FAAH is composed of multiple cavities with mixed hydrophobic and hydrophilic properties, the role of which remains incompletely understood. Anandamide is thought to enter the active site through a "membrane-access" (MA) channel and position its flexible fatty acyl chain in a highly hydrophobic "acyl chain-binding" (AB) cavity to allow for hydrolysis to occur. Using microsecond molecular dynamics (MD) simulations of FAAH embedded in a realistic membrane/water environment, we show now that anandamide may not lock itself into the AB cavity but may rather assume catalytically significant conformations required for hydrolysis by moving its flexible arachidonoyl tail between the MA and AB cavities. This process is regulated by a phenylalanine residue (Phe432) located at the boundary between the two cavities, which may act as a "dynamic paddle." The results identify structural flexibility as a key determinant by which FAAH recognizes its primary lipid substrate.

Characterization of molecular determinants of the conformational stability of macrophage migration inhibitory factor: leucine 46 hydrophobic pocket.

Macrophage Migration Inhibitory Factor (MIF) is a key mediator of inflammatory responses and innate immunity and has been implicated in the pathogenesis of several inflammatory and autoimmune diseases. The oligomerization of MIF, more specifically trimer formation, is essential for its keto-enol tautomerase activity and probably mediates several of its interactions and biological activities, including its binding to its receptor CD74 and activation of certain signaling pathways. Therefore, understanding the molecular factors governing the oligomerization of MIF and the role of quaternary structure in modulating its structural stability and multifunctional properties is crucial for understanding the function of MIF in health and disease. Herein, we describe highly conserved intersubunit interactions involving the hydrophobic packing of the side chain of Leu46 onto the ß-strand ß3 of one monomer within a hydrophobic pocket from the adjacent monomer constituted by residues Arg11, Val14, Phe18, Leu19, Val39, His40, Val41, Val42, and Pro43. To elucidate the structural significance of these intersubunit interactions and their relative contribution to MIF's trimerization, structural stability and catalytic activity, we generated three point mutations where Leu46 was replaced by glycine (L46G), alanine (L46A) and phenylalanine (L46F), and their structural properties, stability, oligomerization state, and catalytic activity were characterized using a battery of biophysical methods and X-ray crystallography. Our findings provide new insights into the role of the Leu46 hydrophobic pocket in stabilizing the conformational state of MIF in solution. Disrupting the Leu46 hydrophobic interaction perturbs the secondary and tertiary structure of the protein but has no effect on its oligomerization state.

Integrating computational methods to retrofit enzymes to synthetic pathways.

Microbial production of desired compounds provides an efficient framework for the development of renewable energy resources. To be competitive to traditional chemistry, one requirement is to utilize the full capacity of the microorganism to produce target compounds with high yields and turnover rates. We use integrated computational methods to generate and quantify the performance of novel biosynthetic routes that contain highly optimized catalysts. Engineering a novel reaction pathway entails addressing feasibility on multiple levels, which involves handling the complexity of large-scale biochemical networks while respecting the critical chemical phenomena at the atomistic scale. To pursue this multi-layer challenge, our strategy merges knowledge-based metabolic engineering methods with computational chemistry methods. By bridging multiple disciplines, we provide an integral computational framework that could accelerate the discovery and implementation of novel biosynthetic production routes. Using this approach, we have identified and optimized a novel biosynthetic route for the production of 3HP from pyruvate.

Pushing the frontiers of first-principles based computer simulations of chemical and biological systems.

The Laboratory of Computational Chemistry and Biochemistry is active in the development and application of first-principles based simulations of complex chemical and biochemical phenomena. Here, we review some of our recent efforts in extending these methods to larger systems, longer time scales and increased accuracies. Their versatility is illustrated with a diverse range of applications, ranging from the determination of the gas phase structure of the cyclic decapeptide gramicidin S, to the study of G protein coupled receptors, the interaction of transition metal based anti-cancer agents with protein targets, the mechanism of action of DNA repair enzymes, the role of metal ions in neurodegenerative diseases and the computational design of dye-sensitized solar cells. Many of these projects are done in collaboration with experimental groups from the Institute of Chemical Sciences and Engineering (ISIC) at the EPFL.

Mixed quantum mechanical/molecular mechanical (QM/MM) simulations of adiabatic and nonadiabatic ultrafast phenomena.

A thorough theoretical description of ultrafast phenomena that occur in complex systems constitutes a formidable challenge. It not only necessitates the use of quantum mechanical methods that can describe ground and possibly even electronically excited state potential energy surfaces with sufficient accuracy but also calls for approaches that can take the real-time dynamics of a system and the coupling between its electronic and nuclear degrees of freedom fully into account. Over the last years, our group has been active in the development of mixed quantum mechanical/molecular mechanical (QM/MM) methods for the in situ simulations of dynamical phenomena in ground and excited states within the adiabatic (Born-Oppenheimer) approximation. Recently, we have extended our theoretical tools with the explicit inclusion of nonadiabatic effects in the framework of Ehrenfest dynamics and Tully's fewest switches surface hopping. These extensions allow the theoretical description of nonadiabatic ultrafast phenomena in the gas phase as well as in solution, and complex biological environments.

Predicting novel binding modes of agonists to ß adrenergic receptors using all-atom molecular dynamics simulations.

Understanding the binding mode of agonists to adrenergic receptors is crucial to enabling improved rational design of new therapeutic agents. However, so far the high conformational flexibility of G protein-coupled receptors has been an obstacle to obtaining structural information on agonist binding at atomic resolution. In this study, we report microsecond classical molecular dynamics simulations of ß(1) and ß(2) adrenergic receptors bound to the full agonist isoprenaline and in their unliganded form. These simulations show a novel agonist binding mode that differs from the one found for antagonists in the crystal structures and from the docking poses reported by in silico docking studies performed on rigid receptors. Internal water molecules contribute to the stabilization of novel interactions between ligand and receptor, both at the interface of helices V and VI with the catechol group of isoprenaline as well as at the interface of helices III and VII with the ethanolamine moiety of the ligand. Despite the fact that the characteristic N-C-C-OH motif is identical in the co-crystallized ligands and in the full agonist isoprenaline, the interaction network between this group and the anchor site formed by Asp(3.32) and Asn(7.39) is substantially different between agonists and inverse agonists/antagonists due to two water molecules that enter the cavity and contribute to the stabilization of a novel network of interactions. These new binding poses, together with observed conformational changes in the extracellular loops, suggest possible determinants of receptor specificity.

Role of aggregation in rhodopsin signal transduction.

Many G protein-coupled receptors (GPCRs) are known to form dimers or even oligomers, and these aggregated states have been proposed as functional units responsible for signal transduction and G protein activation. However, the nature of their involvement has remained elusive. Here, we have investigated the role of aggregation in the signal transduction for dimeric forms of the prototypical GPCR rhodopsin using molecular dynamics simulations. The early steps after photoexcitation are characterized by a tandem mechanism in which one monomer is responsible for light detection while the other serves as the G protein activation site. Dimerization ensures efficient cross-talk between the two units within a few tens of nanoseconds following photoexcitation. This interface-mediated pathway suggests oligomerization-aided signal transduction as a crucial biological strategy to enhance activation efficiency across the entire family of GPCRs.

A conserved protonation-induced switch can trigger "ionic-lock" formation in adrenergic receptors.

The mechanism of signal transduction in G-protein-coupled receptors (GPCRs) is a crucial step in cell signaling. However, the molecular details of this process are still largely undetermined. Carrying out submicrosecond molecular dynamics simulations of beta-adrenergic receptors, we found that cooperation between a number of highly conserved residues is crucial to alter the equilibrium between the active state and the inactive state of diffusible ligand GPCRs. In particular, "ionic-lock" formation in beta-adrenergic receptors is directly correlated with the protonation state of a highly conserved aspartic acid residue [Asp(2.50)] even though the two sites are located more than 20 A away from each other. Internal polar residues, acting as local microswitches, cooperate to propagate the signal from Asp(2.50) to the G-protein interaction site at the helix III-helix VI interface. Evolutionarily conserved differences between opsin and non-opsin GPCRs in the surrounding of Asp(2.50) influence the acidity of this residue and can thus help in rationalizing the differences in constitutive activity of class A GPCRs.

Observation of "ionic lock" formation in molecular dynamics simulations of wild-type beta 1 and beta 2 adrenergic receptors.

G protein coupled receptors (GPCRs) are a large family of integral membrane proteins involved in signal transduction pathways, making them appealing drug targets for a wide spectrum of diseases. The recently crystallized structures of two engineered adrenergic receptors have opened new avenues for the understanding of the molecular mechanisms of action of GPCRs. Taking the two crystal structures as a starting point, we carried out submicrosecond molecular dynamics simulations of wild-type beta(1) and beta(2) adrenergic receptors in a lipid bilayer under physiological conditions. These simulations give access to structural and dynamic properties of the receptors in pseudo in vivo conditions. For both systems the overall fold properties of the transmembrane region as well as the binding pocket remain close to the crystal structure of the engineered systems, thus suggesting that the ligand binding mode is not affected by the introduced modifications. Both simulations indicate the presence of one or two internal water molecules absent in both crystal structures and essential for the stabilization of the binding pocket at the interface between transmembrane helices III, IV, and V. The different interactions arising from the substitution of Tyr308 in beta(2)AR into Phe325 in beta(1)AR induce different conformations of the homologous Asn(6.55) inside the binding pockets of the two receptors, suggesting a possible origin of receptor specificity in agonist binding. The equilibrated structures of both receptors recover all of the previously suggested features of inactive GPCRs including formation of a salt bridge between the cytoplasmatic moieties of helices III and VI ("ionic lock") that is absent in the crystal structures.

Topologically Based Multipolar Reconstruction of Electrostatic Interactions in Multiscale Simulations of Proteins.

We present a new method to incorporate electrostatic interactions in coarse-grained representations of proteins. The model is based on a topologically reconstructed multipolar expansion of the all-atom centers of charge, specifically of the backbone dipoles and the polar or charged side chains. The reliability of the model is checked by studying different test cases, namely protein-cofactor/substrate interactions, protein large conformational changes, and protein-protein complexes. In all cases, the model quantitatively reproduces the all-atom electrostatic field in both a static and a dynamic framework. The model is of general applicability and can be used to improve both full coarse-grained simulations and hybrid all-atom/coarse-grained multiscale approaches.

Microseconds dynamics simulations of the outer-membrane protease T.

Conformational fluctuations of enzymes may play an important role for substrate recognition and/or catalysis, as it has been suggested in the case of the protease enzymatic superfamily. Unfortunately, theoretically addressing this issue is a problem of formidable complexity, as the number of the involved degrees of freedom is enormous: indeed, the biological function of a protein depends, in principle, on all its atoms and on the surrounding water molecules. Here we investigated a membrane protease enzyme, the OmpT from Escherichia coli, by a hybrid molecular mechanics/coarse-grained approach, in which the active site is treated with the GROMOS force field, whereas the protein scaffold is described with a Go-model. The method has been previously tested against results obtained with all-atom simulations. Our results show that the large-scale motions and fluctuations of the electric field in the microsecond timescale may impact on the biological function and suggest that OmpT employs the same catalytic strategy as aspartic proteases. Such a conclusion cannot be drawn within the 10- to 100-ns timescale typical of current molecular dynamics simulations. In addition, our studies provide a structural explanation for the drop in the catalytic activity of two known mutants (S99A and H212A), suggesting that the coarse-grained approach is a fast and reliable tool for providing structure/function relationships for both wild-type OmpT and mutants.